Evolution of Mixed Particles Interacting with Classical Sources 
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We study the systems of scalar and spinor particles with mixing emitted by external classical 
sources. The particles wave functions exactly accounting for external sources are obtained directly 
from the Lorentz invariant wave equations in (3+l)-dimensional space-time. Then we discuss sources 
which are localized in space and emit harmonic radiation. We obtain that the considered scalar 
and spinor fields can be converted from one type to another due to the presence of the vacuum 
mixing. This phenomenon is shown to be analogous to neutrino flavor oscillations in vacuum since 
the calculated transition and survival probabilities coincide with the corresponding expressions for 
neutrino oscillations. When we discuss the evolution of spinor particles, the situation of general mass 
matrix, which has both Dirac and Majorana mass terms, is studied. In this case we analyze the 
appearance of antiparticles in the initial particles beam. The relevance of the developed technique 
to the description of neutrino flavor oscillations is examined. 
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I. INTRODUCTION 

The studying of the evolution of elementary particles 
systems with mixing is an important problem of high 
energy physics. Nowadays we have received the experi- 
mental confirmation that there exists the mixing in the 
lepton sector of the standard model (see, e.g., Ref. [l|) 
which leads to the transitions between different neutrino 
flavors, i.e. to neutrino flavor oscillations. 

Recently considerable efforts were made to theoreti- 
cally describe the evolution of flavor neutrinos with mix- 
ing in order to substantiate the results of the quantum 
mechanical approach to the problem in question (see 
Ref. [HQ). The issue of the description of neutrino flavor 
oscillations with help of field theory methods is important 
since some problems arise when the quantum mechanics 
is used to describe neutrino oscillations. For instance, 
it is not clear (i) whether energies or momenta of dif- 
ferent neutrino eigenstates are equal, (ii) whether neu- 
trino oscillations occur in space or in time, (hi) whether 
it is important to take into account the multicomponent 
structure of neutrino wave functions etc. Some of the 
unresolved questions of the neutrino flavor oscillations 
problem are also briefly outlined in Ref. [J| . 

It is necessary to mention that in Refs. 0, H, 0, H, H, 
liH [Til [T2I [Tl| the quantum field theory was applied for 
the studying of neutrino flavor oscillations in vacuum. In 
Ref. [f| the authors studied neutrino oscillations with the 
most general neutrino mass matrix and derived sum rules 
in neutrino oscillations. The authors of Refs. [6|, [7[ dis- 
cussed neutrino oscillations in a similar way as in Ref. Q . 



'Electronic address: dvmaxim@cc.jyu.fi 



In those works the evolution of neutrino mass eigenstates 
was considered, i.e. one formulates neutrino oscillations 
problem without resort to flavor eigenstates. It should 
be noticed that, in a sense, the problem of neutrino flavor 
oscillations is incorrectly formulated because in various 
processes where a neutrino participates, emitted parti- 
cles can belong to several neutrino generations. For ex- 
ample, together with muon neutrinos in a pion decay, a 
small number of electron neutrinos is produced (see also 
Sec. lIVp . Hence we should rewrite neutrino emission and 
detection vertices in terms of mass eigenstates. Then we 
can examine neutrino emission and detection processes 
by studying incoming and outgoing charged leptons as in 
Refs. 0, H, 0|- These particles are not mixed and thus 
they can serve as properly defined in and out states. 

In Ref. [H the corrections to the quantum mechan- 
ical transition probability formula were obtained when 
the evolution of Fock states of flavor neutrinos was ex- 
amined. In Ref. [9( it was demonstrated that neutrino 
oscillations can be described in terms of entangled flavor 
states. Several processes involving neutrino weak states, 
constructed in Ref. Q , were calculated in Ref. [ToJ • The 
incompatible results were revealed to appear if one uses 
the inequivalent vacua model proposed in Ref. 0. The 
S'-matrix approach was used in Ref. [llj to examine the 
coherence problem in neutrino flavor oscillations. It was 
shown in Ref. [T2| that neutrino oscillations could be sup- 
pressed in a gedanken experiment where a neutrino is 
produced in a crystal by a relativistic electron. Neutrino 
flavor oscillations were analyzed in Ref. fl3j using the 
path integral method. Finally it should be mentioned 
that various theoretical approaches to neutrino flavor os- 
cillations were reviewed in Ref. [T|| which also contains 
the comprehensive bibliography. 

We demonstrated in Ref. [l5j] that neutrino flavor oscil- 
lations in vacuum, in matter and in electromagnetic fields 
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of various configurations could be described in frames of 
the relativistic quantum mechanics approach. In those 
works we formulated the initial condition problem and 
studied the time evolution of neutrino mass eigenstatcs 
using Pauli- Jordan function, or its analog taking into 
account an external field. As a result we received the 
transition probabilities depending on time. Therefore 
that method is the most natural generalization, which ac- 
counts for the coordinate dependence of neutrinos wave 
functions, of the quantum mechanical treatment of neu- 
trino oscillations. The weakness of that approach consists 
in the fact that to obtain a stable oscillations picture one 
should prepare rather broad, in space, initial wave packet 
of a flavor neutrino. Such initial conditions are very diffi- 
cult to implement. In order to overwhelm that difficulty 
one has to consider more realistic models for neutrino 
emission processes. 

This paper continues the series of our publications on 
the theory of neutrino flavor oscillations. The purpose of 
the present work is to provide a deeper understanding of 
neutrino flavor oscillations in vacuum. We study the evo- 
lution of mixed scalar (Sec. [TTJ) and spinor (Sec. IIIip par- 
ticles emitted by classical sources. We start from Lorentz 
invariant Lagrangians in (3 + l)-dimcnsional space-time. 
Therefore all the results are valid for arbitrary energies 
of emitted particles. The fields distributions of scalar 
and spinor particles exactly accounting for the external 
sources are obtained. Then for the external fields with 
the specific time and coordinates dependences we derive 
the expressions for energy densities, for scalar fields, and 
intensities, for spinor fields. These quantities are exam- 
ined in the limit of rapidly oscillating external fields. It 
is shown that the considered expressions coincide with 
the common formulae for transition and survival proba- 
bilities of neutrino flavor oscillations in vacuum. More- 
over in Sec. lIIIBl we study the evolution of fermion fields 
with the general mass matrix and examine the appear- 
ance of antiparticles. We discuss our results in Sec. IIVI 
and suggest that the developed formalism is a theoreti- 
cal model for the process of neutrino flavor oscillations. 
The approximations made in deriving of major results 
are considered in Appendix [A"l 



II. EVOLUTION OF SCALAR PARTICLES 

In this section we study the dynamics of N complex 
scalar particles tp = (tpi, . . . , <{>n) with mixing. The 
scalar particles in question are supposed to interact with 
the external classical fields f\. The Lagrangian for this 
system has the following form: 



A' 



N 



A=l 



AA' = 1 



TV 



E(/W^ a /a), 



(2.1) 



where {m\ x , ) is mass martix. We suppose that it is real 
and symmetric. The non-diagonal elements of this ma- 
trix are called the vacuum mixing terms. The external 
sources f\ = f\(r, t) in Eq. (|2.ip are taken to be arbitrary 
functions of coordinates and time. 

To describe the evolution of the considered system we 
make the matrix transformation of the Lagrangian (|2.1|) 
and introduce the new set of the scalar fields u a , 



N 
a=l 



N 



(2.2) 



in such a way to diagonalize the matrix {m\ x ,), i.e. to 
eliminate the vacuum mixing terms in Eq. (|2.1[) . It 
should be noted that the matrix (U\ a ) in Eq. (12. 2[) is 
orthogonal. 

The Lagrangian expressed in terms of the fields u = 
(til, . • • , un) has the following form: 

N 
a=l 



N 

J2(glu a + ulg a ), 

a=l 



(2.3) 



where m a are the masses of the fields u a . The new sources 
g a are obtained with help of Eq. 



(2.4) 



A=l 



A=l 



The dynamic equation for the field u a is the inhomoge- 
neous Klein-Gordon equation, 



( D - m l) u a = -3a, 



(2.5) 



where □ = —d^d^ is the d'Alembertian operator. It is 
worth mentioning that the squares of masses m 2 , of the 
fields u a are the eigenvalues of the matrix (m\ y ). 

The solution to Eq. (|2.5p for the arbitrary function g a 
has the form (see, e.g., Rcf. 16]) 

u a (r,t) = J d 3 r'dt'D^(v-r',t-t')g a (r',t'), (2-6) 

where D T a et (r,t) is the retarded Green function. This 
function obeys the equation, 



{n-ml)DT{x) = -8 i (x) 1 
D a ct (x) =0 for t < 0, 



(2.7) 



where x* 1 = (t,r). The explicit form of the function 
D a et (x) can be found in Ref. [H, 



D a e \r,t) 



d 4 p 



A=l 



(2tt) 4 m 2 a ~ p 2 + iep° 
i-^){<5( S 2 )-0( S 2 )^J 1 (m aS )}, (2. 
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where 9(t) is the Heaviside step function, Ji(z) is the 
first order Bessel function, s = \Jt 2 — r 2 and r = |r|. 

Now it is necessary to specify the time and spatial co- 
ordinates dependence of the external fields f\(r,t). We 
assume that the sources are localized in space and emit 
harmonic radiation, 



/ A (r,t)=^)/! 0) e MBt <5 3 (r), 



(2.9) 



where /^°' is the amplitude of the source and E is its 
frequency, the analog of the energy of emitted particles. 
In this case Eq. (|2.6[) is rewritten in the following way: 

u a (r, t) = g^e- iEt f drDf (r, r)e iEr , (2.10) 
Jo 



where 



JV 



4 0) = £^/. 



t f(0) 
A ' 



A=l 



is the time independent component of the function g a . 

One can see that two major terms appear in Eq. (|2.10p 
[see also Eq. (|2.8I) ] while integrating over r 



f dTDl ct (r,r)e iET =I 1+ I 2 
Jo 



They are 
1 



2tt 



iEr 9{t - V) j 



Er 



47rr 



(2.11) 



(2.12) 



4?r 



dT6(s 2 ) 



2 Ji(m a s) - Er 



ma ^ r) r dx±t£l J*****, (2.13) 



4-7T 



Vr + 3? 



where x r , 



VP 



and s = yr 2 



It is interesting 



to notice that both X x and T 2 in Eqs. (|2~T2f and (f2TT3|) 
are equal to zero if t < r. It means that the initial 
perturbation from a source placed at the point r = 
reaches a detector placed at the point r only after the 
time t > r. 

Despite the integral I\ in Eq. (|2.12[) is expressed in 
terms of the elementary functions, the integral X 2 in 
Eq. (|2.13|) cannot be computed analytically for arbitrary 
r and t. However some reasonable assumptions can be 
made to simplify the considered expression. An observer 
is at the fixed distance from a source. As we have already 
mentioned one detects a signal starting from t > r. It is 
obvious that a non-stationary rapidly oscillating signal is 
detected when the wave front just arrives to a detector, 
i.e. when t > r. The situation is analogous to waves 
propagating on the water surface. Therefore, if we sup- 
pose that one starts observing particles when the non- 
stationary signal attenuates, i.e. <> c or x m — > oo, we 
can avoid relaxation phenomena. In this case the integral 
T 2 can be computed analytically, 



1 



(e 



iEr 



), 



(2.14) 



where p a = \J E 2 — m 2 is the analog of the particle mo- 
mentum. 

To obtain Eq. (|2.14p we use the known values of the 
integrals, 



J v (mx) 



dx sm(E\/ r 2 + x 2 ) 



\/r 2 + x 

7T 

x J_ 



v/2 



E+y / E T 



dx 



J v {mx) , r-z- — r. 
-ji±=== cos(i?V r 2 + x 2 ) 



\Jr 2 + 



7T 

2^/2 



x N. 



-v/2 



r 



E+yfW 



and the fact that the Bessel and Neumann functions of 
±1/2 order 



J 1/2 (z) =N_ 1/2 (z) 



J-l/2{z) =\l — cosz > 

are expressed in terms of the elementary functions. The 
approximations made in derivation of Eq. (|2.14[) are anal- 
ysed in Appendix [A] 

Using Eqs. (|2~T2]) and (|2~T4|) we obtain the field distri- 
bution of u„ , 



(0) 

^-iEt+ip a r 9a 



f-- (2-15) 
47rr 



Then it is necessary to describe the evolution of the fields 
tp\. With help of Eqs. (|2.2p and (|2.15p one receives the 
field distribution of <p\ in the following form: 



N 



-lEt 



iA'=l 



(0) 
X' 



(2.16) 



We remind that this formula is valid if t ^> r. 

To be more elucidative we consider the system of only 
two scalar fields, (p± and tp 2 . In this case the matrix (U\ a ) 
in Eq. (|2.2p is parametrized with help of the mixing angle 



(Uxa) 



/ COS! 

I sin( 



— sm ( 
cos 8 



(2.17) 



We also adopt the following amplitudes of the external 
sources: f[ ' — and f 2 ^ = f ^ 0. If we discussed 
the evolution of flavor neutrinos, that is in fact done in 
Sec. IIII1 such a choice of the external fields would mean 
that tpi would be a muon or a r-neutrino and (f 2 - an 
electron neutrino. There would be a source of only elec- 
tron neutrinos at r = and a detector at the distance r 
from a source. 
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Using Eqs. (|2.16[) and (|2. 17|) we obtain the fields dis- 
tributions of ipi 2 in the following form: 

pi(r,i) =sin6»cos6»e- ist -^-(e ipir - e ip2r ), 

<p 2 (r, t) = e - iEt -l-(sm 2 9e ipir + cos 2 6e ip2r ). (2.18) 

It was demonstrated in Ref. [l5| that the measurable 
quantity for a scalar field is the energy density rather 
than the field distribution. Constructing the time com- 
ponent T 00 of the energy momentum tensor of <p\ we can 
write down the energy density in the following way: 



W A (r,t) = T > A (r,t)] 

\<P\\~ + |V</?a|~ +m AA |<^ A | 



(2.19) 



Now let us discuss the high frequency approximation of 
the obtained results: E ^> mi. 2. If we again referred 
to the flavor neutrinos evolution (see Sec. IIIip . it would 
mean that ultrarelativistic neutrinos would be emitted. 
Decomposing the "momenta" p a « E — m 2 /(2E) in 
Eqs. (|2.18p and taking into account that H\ « \<p\\ 2 + 
|V<^a| 2 ~ 2E 2 \<px\ 2 one receives the expressions for the 
energy densities 



Hi(r) =2E 



l/l 2 



(47rr) 2 
x i sin 2 (20) sin 2 



(2.20) 



Am 2 
\E ' 



W 3 (r) =2S 2 -]^ 
(47rr) 



x <^ 1 - sin 2 (26») sin 2 



Am 2 
4£ 



where we use the common notation Am 2 



£ 2 



With help of Eqs. (|2.20p we obtain the transition and 
survival probabilities in the following way: 



P. 



alar)/ n . 2(nn\ • 2 f Am 

1 >{r) ~sin (26) sin I — 
f r) (r)~l-sin 2 (2(9)sm 2 ^ 



(2.21) 
(2.22) 



It can be seen that Eqs. (|2.2ip and (|2.22p coincide with 
the usual formulae for transition and survival proba- 
bilities of neutrino flavor oscillations in vacuum that 
were obtained in frames of the quantum mechanical ap- 
proach 0, m. 



III. EVOLUTION OF SPINOR PARTICLES 



The Lorentz invariant vacuum, i.e. without external 
fields, Lagrangian generally has the form, 



A=l 



N 



AA'=1 

+ mfy(vf) T Cvf,+h.c), (3.1) 

where v\' K = (1/2)(1 =p 7 5 )^a are the chiral projections, 
( to aa')> ( m AA')' ( to aa') are * ne Dirac as well as (left 
and right) Majorana mass matrices and C — ij 2 -f is 
the charge conjugation matrix. These mass matrices 
should satisfy certain requirements which are discussed 
in Ref. The Lagrangian (|3.ip is CPT- invariant. In 
this section we adopt the notations for Dirac matrices as 
in Ref. [T^ . The analogous mass terms are generated in 
theoretical models based on the type-II seesaw mecha- 
nism (see, e.g., Ref. [Hj]). 

We supply the Lagrangian (|3.ip with the term describ- 
ing particles interaction with external sources. Suppose 
that a particle is involved in the left-handed currents in- 
teractions, 



E^wi'+^A/r), 



(3.2) 



A=l 



where 7^ = 7 M (1 — 7 5 )/2. The external fields f£ — 
/ A (r,t) are supposed to be arbitrary functions. The 



physical analog of the Lagrangian (|3.2p is the system 
of flavor neutrinos v\ interacting with matter by means 
of the electroweak interactions. In this case t\ is the 
charged SU(2) isodoublet partner of v\. If we study the 
neutrino emission in a process like the inverse /3-decay: 
p + e~ — > u e + n, the external fields in Eq. (|3.2p are (see, 
e.g., Ref. 0) 



f^=-V2G F V n r(l 

fL=o, /& = o, 



a7 5 )^ p , 



(3.3) 



where ^ v and are the wave functions of a proton and 
a neutron, a sa 1.25 and Gf is the Fermi constant. In 
Eq. (|3.3p we assume that a source consists of electrons, 
protons and neutrons. 

After the diagonalization of the whole mass term of 
the Lagrangian (|3.ip one obtains 2N Majorana particles 
with different masses [3j . In the following we discuss the 
case when only Dirac mass term is presented. Then we 
study the general situation. 



In this section we consider the problem analogous to 
that we discuss in Sec. HH namely the evolution of N 
spinor particles u = (v\, . . . , ujsr) with mixing under the 
influence of the external classical fields / A . 



A. Dirac mass term 

In this section we suppose that only Dirac mass matrix 
is non-zero in Eq. (|3.ip . Analogously to the discussion 
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in Sec. [TT] we introduce the new set of spinor fields ip — 
. . . , iPn) by means of the matrix transformation. In 
our situation the mixing matrix is the N x N unitary 
matrix (U\ a ), 

N N 

vx=J2 V****' ^ = Y1 $° U lx- ( 3 - 4 ) 



a=l 



a=l 



The corresponding mass eigenstates ip a are Dirac parti- 
cles. 

When we transform the sum of the Lagrangians (|3.1[) 
and (|3.2[) using Eq. (|3.4|) . it is rewritten in the following 
way: 

N 

£(V0 = ^2^a(i^d^ - m a )ip a 

a=l 

N 

+ + ^**»)> (3.5) 



where \? a is the external source for the fermion tp a , 

N N 
*a = E C/ a t A7^ 1 *a = £^! MC/ Aa. (3.6) 



Using Eqs. (|3.5p and (|3.6p we receive the inhomogeneous 
Dirac equation for the fermion ifj a , 



(i7 M d M - m a )i) a = -* . 



(3.7) 



We mention that the masses m a are the eigenvalues of 
the matrix (m AA ,). 

The solution to Eq. (|3.6[) for the arbitrary spinor \l/ a 
is expressed with help of the retarded Green function for 
a spinor field (see again Ref. [HI]). 

Mr,t)= j 'dVd^ ct (r-r',t-i')* Q (r',i'). (3.8) 

The explicit form of <S£ ct (r, t) can be also found in 
Ref. [iff, 

S^(t,t) = {rfd tl + m a )D^{t,t), (3.9) 

where the function D^ Gt (r, i) is given in Sec. HT1 

To proceed in further calculations it is necessary to 
define the behavior of the external sources. Let us sup- 
pose that the external fields depend on time and spatial 
coordinates in the similar way as in Sec. UH namely 

e x (r, t)/£(r, t) = 6(t)l x f^e- iEt 6 3 (r), (3.10) 

where l\ = is the time independent component of 
the spinor i\ and is the amplitude of the function 

/£(r,i). Using Eqs. (HU) and (f3TT0|) we obtain for * Q , 



N 



(3.11) 



With help of Eqs. (|3~9|) - (l3~TT|) we can rewrite Eq. (f3T8|> 
in the following way: 



d 



^ Q (r,i) = (i 7 ^ — +m a 



x < e 



-\Et 



J*drD^(v,r)e iB ^^\ (3.12) 



Now it is possible to use the technique developed in 
Sec. |TTJ Finally we obtain the field distribution of the 
fermion ip a , 



(3.13) 



where O a = 7°-E — ( 7 n)p a + m a and n is the unit vector 
towards a detector. It should be noted that in deriving 
Eq. (|3.13[) we differentiate only exponential rather than 
the factor 1/r because the derivative of 1/r is propor- 
tional to 1/r 2 . Such a term is negligible at large distances 
from a source. We also remind that Eq. (|3.13|) is valid 
for t ^ r. 

Let us turn to the description of the evolution of the 
fields v\. Using Eqs. (|3.4|) and (|3.13|) . we obtain the 
corresponding field distribution, 



N 



iEt \ * TT. Tit Jp a rrt „,L/. . J\' 



(0)M 



a\' = l 



Airr 



(3.14) 



Analogously to Sec. [TT] we study the evolution of only 
two fermions. The mixing matrix is given in Eq. (|2.17|) . 
We also choose the amplitudes of the sources in the fol- 
lowing way: /f^ = 0, / 2 (0)m = p ^ and l 2 = I. The 
reason for such a choice is similar to that in Sec. HTJ If wc 
suppose that v\ corresponds to a muon or a r-neutrino 
and vi - to an electron neutrino, then it means that there 
is a source of only electron neutrinos and we observe par- 
ticles at the distance r from a source. Using Eqs. (|2.17|) 
and (|3.14p we obtain the fields distributions of each of 
the particles v\ 2 , 



V\ (r, t) = sin 9 cos ( 



AEtF_ 



x=i 



x (e'^Oi - e ip2r 2 )rfl, 

x (sin 2 9e ipir O! + cos 2 Oe^ '0 2 )rfl. (3.15) 

In the following we discuss the high frequency approxi- 
mation: E 3> mi, 2. It again corresponds to the emission 
of ultrarelativistic neutrinos. 

It was shown in our previous works [l5| that the mea- 
surable quantity for a spinor particle is the intensity 
which is defined by the following expression: I\(r,t) = 
\v X (r, t)\ 2 . With help of Eqs. (|3.15p we get the intensities 
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in the form 



(3.16) 



I 2 (r)=-2E 



2 f^r 

{Airr) 2 



Ft*) 



xil- sin 2 (26»)sin 2 



Am 2 
4£ 



r +0 



E 2 



where (T M „) = i^t (0:11)7,, J and a = 7°7. It is possible 
to calculate the components of the tensor (T M „). They 
depend on the properties of the fermion £2 = £, 

(Too) = -H, (Toi) = (T M )* = m~ i[n x Qi, 
(Tij) = %(vn) - (w^nj + Vj^) + ie ijk n k , (3.17) 

where v = (a) is the velocity of fermion £, £ = (S) is its 
spin and S = 7 s a. It should be noted that, if we study 
the evolution of flavor neutrinos and z/2 corresponds to 
an electron neutrino, the fermion t is an electron. 

Let us discuss the simplified case when the spatial com- 
ponents of the four-vector / M are equal to zero: f = 0. 
It corresponds to the neutrino emission by a non-moving 
and unpolarized source. Using Eq. (|3 . 1 7[) we can rewrite 
Eq. (|3.16p in the following way: 

a,_^ l/T „ s _ 3 /o^„ 4 2 Mm 2 



h(r) «2£ 2 (vn) 



h{r) «2£ 2 (vn) 



(47rr) 2 

l/T 



sin 2 (20) sin 2 



4E 



(47rr) 2 
1 - sin 2 (20) sin 2 



4E 



(3.18) 



Now we are able to introduce the transition and survival 
probabilities on the basis of Eq. (|3 . 1 8[) , 



P 2 ( !!i nor) (r)~ sin 2 (20) sin 2 



4E 



„(spinor) / , . 2tnn\ ■ 2 

™2->2 v) ~1 — sm sin 



4£ 



(3.19) 
(3.20) 



It should be noted that Eqs. (|3~TU|) and (j3T2TJ)) are the 
same as the common formulae for the description of neu- 
trino flavor oscillations in vacuum 0, Q . These expres- 
sions also coincide with Eqs. (|2.21[) and (|2.22[) derived for 
the case of the mixed scalar fields. 



B. General mass term 



To begin analyzing the dynamics of mixed fermions 
with the general mass term (|3.ip we rewrite the interac- 
tion Lagrangian (|3.2p in the equivalent form, 



N 



L All , 



A=l 



(3.21) 



Now we can introduce mass eigenstates ip a . The analog 
of Eq. is (see also Ref. 0) 



27V 



o=l 



IN 
o=l 



(3.22) 



where the matrix (U\ a ) is the rectangular N x 2N matrix. 
Note that it is not unitary. We also mention that the 
mass eigenstates ip a are Majorana particles. To imple- 
ment the complete diagonalization of the mass term (|3.ip 
it is necessary to use 2N x 2N unitary matrix, i.e. along 
with Eq. (|3.22[) one should consider also the transforma- 
tion of the right-handed chiral components of the parti- 
cles V\, 



2N 

E^e 

a=l 



(3.23) 



Note that the matrices (U\ a ) and (V\ a ) in Eqs. (|3~22|) 
and (|3.23[) are independent. The modern parameteriza- 
tion for these matrices is given in Ref. [2l| . 

The Lagrangian for the particles ip a takes the form, 



2A 



a=l 



£W)=J2^d^ 

=i 

2N 

E (»"„$W " W^a + h.c.) , (3.24) 



where the external sources ^ a = have the same form 
as in Eq. (|3~6|) . 

Note that Majorana spinors are equivalent to two- 
component Weil spinors [231 ] - Hence we can rewrite the 
spinors ?A^ R and ^f a as 



where 



tit 



( W2VZ 



A 



ba = E Ulxxxfl 



(3.25) 



(3.26) 



A=l 



To obtain Eq. Q3.26P we suppose as in Sec. IIII Al that the 
vectors f£ have only time component /°. We also assume 

that (£\) T = (0,xa)- To derive Eq. (|3~26| it is crucial 
that only left-handed currents interactions are presented 

in Eqs. and (|3l>]) . 

It is useful to rewrite the Lagrangian (|3.24[) in terms 



of the two-component spinors rj = 

(Vl,---,V2N) and (01, ■ • • ,<^2JV) U2 



{'h 



,V2N), V = 



2A 



C( V ,r,*)=J2iriUdt-<TV)r, a 



1 A 



E ( ^ m ^2V* a + vUa + h.c. ) . (3.27) 
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Using Eq. (]3 . 27[) we can receive the wave equation for 
two-component spinors, 



- erV { ;/„ - '/',,T2'/, ! , : : i' 



(3.28) 



The solutions to Eq. (|3.28[) have the form (see, e.g., 
Ref. H), 

Vu(r,t)= J ^v'At'STiv-T'^t-t'^aiv'X), (3.29) 
»£(r, t) = [ d 3 r'd*X ct (r ~r',t- t')^' , f ), (3.30) 



where the retarded Green functions are expressed as (see 
also Ref. HI), 



Kl et (r 7 t) =im a a 2 D^(r,t). 



(3.31) 



Here = (d t , V), = (1, a), 1 is the 2x2 unit matrix 
and DJ t ct (r,t) is given in Eq. (|2.8|) . One can check that 
Eq. (|3~2T)|) and (j330"]) . along with the definition of the 
retarded Green functions (|3.3ip . represent the solutions 
to Eq. (|3.28p by means of direct substituting and using 
Eq. (|2~7D . 

Let us assume that the sources 4>a depend on time and 
spatial coordinates as in previous sections, 



a (r,i) =0(t)^e- [Et d 3 (r), 



N 



(3.32) 



A=l 



where X\ is the time independent component of xa- 
In deriving of Eq. (|3.32p from Eq. (|3.26p we introduce 

the new quantities f\ = f^°, where is the time 

independent part of /°, to simplify the notations. 

On the basis of Eqs. p.29p ~ (|3.32p and using the tech- 
nique of the previous sections we get the particles wave 
functions, 

,(0) 

Va (r,t)=i^d^e- iEt+i ^ r }^-, (3.33) 



,(o) 



(3.34) 



To derive Eqs. (|3.33|) and (|3.34|) we suppose that t^>r. 
The derivatives in Eq. (|3.33p are applied on the exponent 
only because of the same reasons as in Sec. IIII Al 

Using Eqs. l[3~22| . (|3~33| and (j3~34]) as well as the fol- 
lowing identity: 



2 A 



a=l 



(^) C =(^) R = E C/ Aae 



(3.35) 



we get the wave functions of v\ and {y\) as 

2N N 



a=l A' = l 



2N N 

2h : P( \ - x - m„ 



(3.36) 



(«*)" ( r ' ') - - S^ lSt E E _g^ r <Mlv 

a=l A' = l 



where xl? 



(1/2)[1 - (<rn)]x5?. 
Eq. (|3.36l) we suppose that 



(3.37) 



In deriving of 



that is valid for relativistic neutrinos. 

It is possible to construct two four-component Majo- 
rana spinors from two-component Weil spinors, 



Va 



(2) 



( iV2Va 



(3.38) 



where rj a and rj* are defined in Eqs. (|3.33p and (|3.34p . 
One can see that these spinors satisfy the Majorana con- 



dition 



(1,2) 



= tpa ■ We use the spinor ipa' to re- 

_______ ______ 

ceive Eq. (|3.36p and - Eq. (|3.37p . In this case we 
get that v\ and \y\) are obtained as a result of the 
evolution of particles emitted from the same source. 

It should be mentioned that both v\ and {yYf m 
Eqs. (|3.36p and (|3.37p propagate forward in time. To 
explain this fact let us discuss the complex conjugated 
equation (|3.29p . Performing the same computations one 
arrives to the analog of Eq. (|3.37p which would depend 
on time as e lEt . However this wave function describes 
a particle emitted by the source different from that dis- 
cussed here. Indeed, if we studied the complex conju- 
gated Eq. (|3.29p . the integrand there would be S T a et *<j>l. 
It would mean that the source in Eq. (|3.32p would be 

proportional to (f>a e lEt , that, in its turn, would sig- 
nify that iyY) would be emitted in a process involving a 
lcpton which is a charge conjugated counterpart to that 
discussed here. 

Let us illustrate this problem on more physical exam- 
ple. Suppose that we study a neutrino emission in a 
process like the inverse /3-decay, 



2N 

E^a' 
a=l 



(3.39) 



where p, n and £ x stand for a proton, a neutron and for 
a negatively charged lepton. The complex conjugated 
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Eq. (|3.29p would correspond to a process, 

27V 

o=l 
2N 

^P + Y, U *^a, (3-40) 



where v x and l\ denote an antineutrino and a positively 
charged lepton, which is the charge conjugated counter- 
part to £ x . In Eq. (|3.40|) we use the facts that only 
left-handed interactions exist in nature and the fields ip a 
describe Majorana particles. As one can see, Eqs. (|3.39|) 
and (|3.40jl represent two different processes. Hence, if 
we used (r] a )* to obtain {yY) ( r J^); i- e - after a beam of 
neutrinos passes some distance r, it would correspond to 
the initial reaction (|3.40[) rather than (|3.39p . Note that 
the same result also follows from Eq. (|3.38[) if we replace 

Wa <-* ipa there. 

Let us suppose for simplicity that the momentum of 
the fermion £ x is parallel to the neutrino momentum. It 
takes place if a relativistic incoming lepton is studied. 
We also assume that this fermion is in a state with the 
definite helicity, 



i[l-(<rn)] X f =xi 0) . 



(3.41) 



This expression is again natural for the relativistic 
fermion i\. One can notice that interaction La- 
grangian ()3.21|) is written in terms of the left-handed 
chiral projections of l x . Therefore, if we study a rela- 
tivistic lepton, it will have its spin directed oppositely to 
the particle momentum as one can see from Eq. (|3.41[) . 
Using Eqs. ([335)1 . (j33?j) and (|3T4"T|) as well as 

the orthonormality of the two-component spinors X x , 

(x x °^ X^yJ = we get the probabilities to detect v\ 

and (1 



*%) as 

2JV N 

p ^w~E E ei(Pa_Pb)r 

ab=l X' = l 

*u Xa ui x ,ui xl u x , b \h,\ 2 , 

2N N 
ab=l A'=l 

x U* Xa Ul xl U^ xl U x , h \f xl \ 2 



m °- mb r i(Pa-p b )r 

{2Ef 



(3.42) 



(3.43) 



In Eqs. ([3~i2]) and ([3~l3"]) we drop the factor 

(2.E)7(47rr) 2 . 

We can see that Eqs. (|3.42j) and (|3.43|) contain the os- 
cillating exponent. This our result reproduces the usual 
formulae for neutrino oscillations in vacuum Q . It should 
be also noticed that the expressions for the probabilities 
depend on r rather than on t in contrast to our previous 
works 15]. Note that the problem whether neutrino os- 
cillations happen in space or in time was also discussed 



in Ref. [24J . The similar coordinate dependence of the 
probabilities was obtained in fiefs. 0, ||| @] where the 
problem of neutrino oscillations in vacuum was studied. 

It was mentioned in Ref. [5] that oscillations between 
active and sterile neutrinos are possible in case of the 
non-unitary matrix (U Xa ), i.e. the presence of only Ma- 
jorana mass terms is not sufficient for the existence of 
this kind of transitions. The situation is analogous to 
that considered in the pioneering work [2{| where oscilla- 
tions between neutrinos and antineutrinos were studied. 
In Eq. (|3.43p we obtain that the probability to detect 
(1%) is suppressed by the factor m a rrn,/E 2 . It was also 
mentioned in Rcfs. [l,[26j]. 



IV. DISCUSSION AND CONCLUSIONS 

In conclusion we mention that we studied the evolu- 
tion of the systems of mixed particles emitted by clas- 
sical sources. The cases of both scalar and spinor par- 
ticles were considered. The fields distributions of parti- 
cles were found for the arbitrary sources dependences on 
time and spatial coordinates. Note that the expressions 
for the fields distributions were received directly from the 
Lorentz invariant wave equations and thus they are valid 
for arbitrary energies of emitted particles. It should be 
also mentioned that all the calculations were carried out 
in (3 + 1) dimensions. 

We demonstrated that the expressions for energy den- 
sities, in the scalar fields case, and intensities, in the 
spinor fields case, coincide with the usual transition 
and survival probabilities of neutrino flavor oscillations 
in vacuum derived within the quantum mechanical ap- 
proach. Since we study the process of neutrino oscilla- 
tions on the basis of the field theory method we repro- 
duce not only the quantum mechanical probabilities as 
the leading terms, but also one can derive the correc- 
tions to these expressions. It is known (see Refs. @,[ll) 
that the Pontecorvo formula for neutrino oscillations has 
the corrections ~ (m a /E) 2 , that are small for ultrarel- 
ativistic neutrinos. In the present work we showed that 
the analogous terms also appear in the expressions for 
probabilities [see, e.g., Eqs. (|2.20[) and (|3.16p ]. In ad- 
dition, in Sec. IIII Bl we considered the possibility of the 
appearance of antiparticles in the initial particles beam. 

In our previous works 



151 ] we studied neutrino fla- 



vor oscillations in vacuum and in various external fields 
within the field theoretical approach. In those papers we 
solved the initial condition problem for the wave equa- 
tions describing the evolution of neutrino mass eigen- 
states. To study the dynamics of the particles one should 
use Pauli- Jordan function as the Green function for the 
considered initial condition problem. As a result we re- 
ceived the final wave functions corresponding to the time 
t of the system evolution. It was shown that the corre- 
sponding probabilities to detect a flavor neutrino coin- 
cide with usual probabilities of neutrino oscillations de- 
rived in the Schrodinger equation approach. Therefore 
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we obtained that neutrino oscillations happen in time. 
The weakness of that method was that to get the stable 
oscillations picture, i.e. reproduce the Pontecorvo for- 
mula, one should prepare very broad initial wave packet. 
Physically it would correspond to a permanent, or rather 
long-continued, neutrino source. 

One can also treat neutrino flavor oscillations in vac- 
uum describing the propagation of neutrino mass eigen- 
states as the internal lines in a Feynman diagram [!, [a 0] • 
Thus a Feynman propagator serves as a Green function in 
that approach. In that method one replaces flavor eigen- 
states with mass eigenstates at the neutrino emission and 
detection points. It is possible to say that one describes 
neutrino oscillations without resort to flavor neutrinos. 
It was found that transition and survival probabilities 
depend on the distance between emission and detection 
points. Although that approach seems to be quite natural 
it encounters some problems. For example, it is unclear 
(i) if different mass eigenstates should have equal ener- 
gies, (ii) if one loses coherence when the energy is mea- 
sured with high accuracy, (iii) if the usual Pontecorvo 
formula is unique because the non-standard oscillations 
formulae are possible in that approach (see the criticism 
in Ref. [3). 

In the present work we developed another approach to 
the neutrino flavor oscillations problem. In frames of our 
method we derived the coordinate dependent transition 
and survival probabilities. We examined the relaxation 
phenomena occurring when a neutrino signal arrives to 
a detector. It was revealed that the stable oscillations 
picture could be implemented if the time after the neu- 
trino emission is much greater than the distance between 
a source and a detector. It means that the quantum 
mechanical transition probability formula is valid for a 
rather long-continued neutrino source. This result, in a 
sense, confirms the results of our previous works [15f (see 
also above). Note that one can come the similar conclu- 
sions basing on Ref. [13] where neutrino oscillations were 
studied on the basis of an exactly solvable model. 

The weakness of the present method is that in various 
processes where a neutrino participates, particles belong- 
ing to different flavors can be emitted (see, e.g., Ref. [28f). 
For instance, there are the following decay channels of a 
charged pion: n + — > /i + + and 7r + — > e + + v e . There- 
fore in frames of our model the classical sources of both 
muon and electron neutrinos would correspond to this 
process. However the latter channel is suppressed by the 
factor 10 4 [29]. Thus we can suppose that only one type 
of neutrino sources is presented. Nevertheless one can 
formally consider the sources of all neutrino species in 
frames our method. The amplitudes of the sources 
should be taken either from the quantum field theory of 
the neutrino emission or from an experiment. 

Note that neutrino flavor oscillations in the model with 
a source and a detector were described in Ref. [3(J. In 
that paper the source of particles was also assumed to 
be spatially localized. However its time dependence was 
chosen to be gaussian that is different from our approach. 



Summarizing we can say that our method for the de- 
scription of the flavor neutrinos evolution can be useful 
when we have spatially localized neutrino sources which 
start to emit neutrinos at some moment of time. One 
of the possible examples is the neutrinos emission in a 
supernova explosion. 
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APPENDIX A: EVALUATION OF INTEGRALS 

It is interesting to evaluate the inexactitude which 
is made when we approach to the limit x m — > oo in 
Eq. (|2.13p . Let us discuss two functions 



F(r,t) 



' dx J l( mX ) c iEV^+^ 
ii \Jr 2 + x 2 

"d^^, (Al) 



P 



(A2) 



These functions are proportional to I2 in Eqs. (|2.13|) 
and (|2.14p respectively. In Eqs. (|A1|) and (|A2jl we use 
dimensionless parameters p = mr, £ — Er = 7/5, 
7 = E/m, y m = V{t/r) 2 - 1 and P = y/£ 2 -p>. 

In case we study neutrinos, we get that £ ^> p » 1 
in almost all realistic situations. For example, suppose 
we study a neutrino emitted in a supernova explosion is 
our Galaxy. The typical distance is r ~ lOkpc. Taking 
to ~ 1 eV and E ~ 10 MeV, we receive that £ ~ 10 35 
and p ~ 10 28 . 

Basing on the analysis of Sec. |TT] we can rewrite 
Eq. ([AT]) as 



F(r,t) =F (r)-6F, 

p + OC 

5F= I dy 



VTTy 2 



(A3) 



Using the fact that p ^> 1, y m 3> 1 and the representation 
for the Bessel function, 



JAz) ps \ — cos ( z — — ) at z — > +00, (A4) 

V 7TZ \ 4 ' 



we obtain for the function SF in Eq. (IA3|) the following 
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FIG. 1: The absolute values of the functions F(r,i) and Fo(r) 
versus t. 



expression: 
SF «- 



= {ci([7 + l]py m ) + ci([ 7 - l]py m ) 



+ iN([7 + l]pWm) + si([7 - l]Wm)]}, (A5) 



where ci(z) and si(z) are cosine and sine integrals. Using 
the asymptotic expression, 



|ci(z)| ~ |si(*)| ~ - 



at 



+00, 



(A6) 



we obtain that the function 5F approaches to zero as 
l/(y m p 3 / 2 ) at great values of y m and p. Note that this 
result remains valid for a particle with an arbitrary 7 fac- 
tor, i.e. rapid oscillations of the function 5F will atten- 
uate even for slow particles. This analysis substantiates 
the approximations made in Sec. HU 

Finally let us illustrate the behavior of the functions 
F(r,t) and Fo(r). On Fig. Q] we present the absolute 
values of these functions versus t. This figure is plotted 
for p = 1 and £ = 10. The solid line is the absolute value 
of the function F(r,t) and the dashed line - F (r). As 
we mention in Sec. [TT] the relaxation phenomena occur 
when t > r. It can be also seen on Fig. [TJ It is possible 
to notice that \F(r, t) \ — > l-Folj)! at great values of t as it 
is predicted in Sec. [Til 
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